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Abstract 



We study transport in a class of exactly solvable models of interacting fermions in one di- 



O ' mension. We contrast these models with models of non-interacting fermions in an Aharanov- 

■ Bohm ring to which they are superficially similar. We introduce magnetic and non-magnetic 

impurities at a site, through either a weak ^-function potential or through a weak link. Using 

^ ! a renormalisation group analysis, we show that the strength of the nonmagnetic impurity is 

O \ not affected by the interaction, whereas the magnetic impurity cuts the wire at the impurity 

,__! ■ site. 
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Strongly correlated transport is an emerging field, which has suddenly shot into promi- 
nence in the last few years. One reason has been the advances in semiconductor technology, 
which have enabled the fabrication of one- dimensional quantum wires in the single channel 
limit |l[]. Besides, the recent discovery of carbon nanotubes0 and in particular, the actual 
fabrication of single- wall carbon nanotubes and multi probe experiments [§] on them have 
spurred both theoretical and experimental activity in this field. 

One dimensional interacting electrons are known to exhibit Luttinger liquid (LL) be- 
haviour, rather than Fermi liquid behaviour, characterised by spin-charge separation and 
absence of Landau quasi-particles[|]. In the presence of impurities, transport responses are 
marked by interaction- dependent power laws. The dramatic effect of interactions on trans- 
port in a LL was shown in a paper by Kane and Fisher]^, who studied transport in LLs, 
through both single and double barrier structures. They showed through renormalisation 
group arguments that spinless electrons with repulsive interactions, incident upon a single 
barrier are completely reflected at zero temperature and for vanishing bias, whereas for at- 
tractive interactions there is perfect transmission. This analysis has been further modified 
by several groups P] who included other corrections, but essentially verified their analysis. It 
has also been verified by an exact solution at arbitrary coupling constant using the thermo- 
dynamic Bethe ansatz0. Furthermore, their analysis has also been extended to other cases, 
such as with addition of external biases P], crossed LLsP] and double and multiply crossed 
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All of these analyses, are essentially within the bosonised Luttinger model approach, with 
strong or weak barriers, whose effects are then extrapolated using renormalisation group 
arguments. Since the difference between transport for the non-interacting Fermi liquid and 
the interacting LL is so dramatic, it is worthwhile to see whether such results are reproduced 
in other models of interacting fermions. With a view to addressing this problem, in this 



paper, we study transport in a specific exactly solvable model |]TT| of interacting fermions in 
one dimension. This model is described in terms of two species of fermions, with pseudospin 
index a = ± and with 'gauge' interactions described by a Hamiltonian given by 



N 



^ = EE«/(n..r- (i: 



Here, IIo-j = Po-i + cr^o-(a^o-i) is the 'covariant momentum' introduced in Ref.[|Tl[ and A^ is a 



'band index' introduced in Ref. |T^ r^|, which generalised these models by allowing for higher 
powers of the covariant momenta in the Hamiltonian. Particles interact with each other via 
the gauge potential, given for the particle at the position x by A„{x) = J2j V{x — X-^j) 
- i.e., the potential for a particle with positive pseudospin is due to the presence of the 
particles with negative pseudospin and vice-versa. Note that the potential depends on all 
the particles of the opposite sign of pseudospin irrespective of the number of bands. The 
potential is chosen to be an even function, vanishes at infinity and explicitly breaks time- 
reversal invariance, although it is invariant under a combined operation of time reversal and 
reversal of pseudo-spin index. In this paper, we shall confine ourselves to the original model 
of Schulz and Shastry|Tll| and set A^ = 1 and 1 = 1. 



To study the problem of transport in this model, we first note that we cannot study the 
model on an open wire as is normally done for transport problems, because here interactions 
are essentially introduced through a magnetic flux. Particles of one species (or pseudospin 
orientation) give rise to an effective Aharanov-Bohm flux acting on the particles of the other 
species(other pseudospin orientation). Such interactions cannot be introduced in an open 
wire. Hence, to study equilibrium transport in these models, we shall consider fermions 
on a ring of length L -i.e., we impose the periodic boundary conditions iIj{x) = ip{x + L). 
Instead of measuring current in the limit of vanishing external voltage, we shall measure the 
persistent current in the model in the limit of vanishing external flux through the ring. 

By performing a pseudo-unitary transformation on the model, 

g^5({.+.},{-.})p^^e-^^«^+^>-^^-» = P.^ - 9..,5({x+,}, {x,}), (2) 

the interaction term in Eq.(|l|), is eliminated and the Hamiltonian is transformed to the free 
one given by if = (1/2) ^aiiVcriY ■ However, the boundary conditions on the wave-functions 
are now different and the quantisation conditions for the momenta are now replaced by 

Lk^i ip N^5 = 27in±^i. (3) 

Here 5 is a phase shift that can be computed in terms of the unitary transformation func- 
tion S, which in turn is related to the original interaction, since S is chosen to cancel the 



interaction |]TT|. n±^i are quantum numbers analogous to those used in the non-interacting 
case. Note that 6, which is defined as the fractional part of ^A^± (5 takes value from 
—.5 to +.5) leads to interactions; the integer part of the ostensible interaction term ^N± 
merely gives integer changes in the quantum numbers n±^i, which anyway take values over 
all integers. 

When there are no interactions, - i.e., when we have free electrons on the ring(or ^A^± is 
an integer), - and when there is no external flux through the ring, clearly there is no current. 
However, once we introduce interactions, either by changing 6 or the number of particles, 
essentially, an intrinsic flux is introduced through the ring, which causes a bias. Thus, if we 
plot the persistent current as a function of the external flux as is normally done, we find 
a shift in the plot with the maximum of the current occurring when 0ext = n(f)Q/2 — 6 and 
the minimum or zero current occurring at (pext = n(f)Q — 6, where n = integer and 0o is the 
flux quantum. This implies a non-zero magnetic moment of the interacting ring, even in the 
absence of any external driving flux. In this model, this is not surprising, since the model 
intrinsically violates the discrete time reversal (TR) symmetry. 

Our aim is to study transport in this theory at low temperatures and at low external 
driving forces (here an external flux), when an impurity is introduced in the model. Since for 
open wires, interactions dramatically change the results we expect the differences between 
interacting and non-interacting fermions on the ring also to be dramatic for transport through 
barriers or constrictions. How do we see this ? Once any impurity is introduced at a site, 
the model ceases to be exactly solvable. Since one has a first quantised formulation of the 
model, one may expect that a Landauer-type scattering matrix approach to transport (which 



is applicable for free fermions) may be a viable way to study impurity scattering, although 
here the fermions are interacting. However, we find that the infrared catastrophe, which 
changes the states of all the particles in the system whenever a single particle is added to or 
removed from the system, prevents a convenient analysis in the first quantised formulation. 
In other words, because of the strong interaction, which implies a collective motion of all the 
particles, it is impossible to isolate a single particle, study its transmission and reflection, 
and then add up the contributions of all the particles as is done in the Landauer-Buttiker 



formalism. Hence, we follow the more standard procedure of first bosonising[|Tj] the low 
energy (collective) excitations, then rewriting the impurity scatterer in terms of the bosonic 
variables and studying its effect perturbatively, and finally using renormalisation group to 
extrapolate the results. 

We first obtain the ground state energy Eq in a sector with N± particles by choosing n±^i = 
^±i =*= [[^^T]]int7 where n!^- are the quantum numbers in the absence of any interaction. 
Excitations about the ground state are obtained by constructing the low energy Hamiltonian 
when particles are added to the left and right Fermi points. Let us assume that the ground 
state has 2no + 1 particles of each kind and that ^{2no + 1) is an integer. The addition of 
n±R{n±L) particles at the right(left) Fermi points will cause a second order change in the 
ground state energy given by 
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where the initial particle density is denoted as n = ^ "^"^ , 6 = ^, and N± = n±R + n±L 
and J± = n±ji — n±L are the particle and current quantum numbers respectively. 

Let us introduce boson fields 0± and their conjugate momenta H-t = c?t-0± = dx9±, 



following the notation of Shankar [^ . These are related to the current and charge densities 
as 

N± = —=d^(t)± and J± = ^9^6*^. (5) 

'TX \ TX 



So, the effective low energy Hamiltonian including quantum fluctuations is now given by 

H = —jdx{il + 5')[{d^<P+f + d,<pl)] + [{dM' + idj.f] + 25idJ^d.,<P^ - dJ+d,<P^)}. 

(6) 



e+±e. 



In terms of the standard spin and charge fields defined by 9c/ s = '~fK~ ^^'^ *i^c/s = /t ■• 
the Hamiltonian can be rewritten as 

777r C — — 

H = —J dx{{l + 6^)[{d,^cf + d,(Pl)] + [{dJcf + {dJs?] + 26{dJA(Ps - dJsd.cPc)}. (7) 

Note that unlike the usual models of LLs, here, the Hamiltonian is not separable in terms 
of the spin and charge variables - the S term couples them, just as it couples the ± fields. 
This is yet another motivation to study transport in this model. Until now, the Kane-Fisher 
results]^ have only been obtained in models where spin and charge are explicitly separable. 
Note also that we may work either with spin and charge fields {a = c/s in Eg . (p!0D ) or with 
the + and — fields (a = ± in Eq.(|l^)). Both representations are equivalent. 



The Hamiltonian can be brought to a form where it looks almost non-interacting by 
defining new fields 

0± = dj± =F 5(f)^ and 0± = 0±, (8) 

or Oc/s = d^dc/s ± S(f)s/c and 0c/s = 4>c/s- (9) 

The new fields 6± and (f)± {9c/s and 0c/s) satisfy the same commutation relations as the 
non-tilde fields, - [6'q,(x), 0^(1/)] = iQ{x — y) = [6'q,(x), 0q(i/)]. In terms of these fields, the 
Hamiltonian apparently looks spin-charge (or equivalently +/— ) separable and is given by 

^ = X ^ jdx{{d^~^^f + {dj^f}. (10) 

Q=it or c/s 

Either one of the fields can be treated as the coordinate and the gradient of the other field 
can be treated as the momentum. 

We now study the effects of impurity scattering in this bosonised model. We will work 
in the limits of a weak barrier and a weak link and then extrapolate to intermediate values 
of the coupling strength of the impurity. 

Let us first consider scattering from a small barrier. A potential scatterer or impurity in 
the model can be introduced as 

<5iJ = \/(V'l(0)^+(0)+ ^1(0)^40)) (11) 

where the first and second terms backscatter + and — particles respectively. In momentum 
space, this term translates into 2kp (and higher momentum transfer) scatterings between the 
left and right Fermi points for both the + and — particles. The lowest order 2kF scattering 
term is given by 

5H^ = VeY. ^l,Li^^,R + h-c. (12) 

where either a + or — particle is back scattered. Here, both charge and spin degrees of 
freedom are involved. However, we can also consider higher order scatterings {e.g., Akp 
scattering) where both + and — particles are scattered such as 

6H2 = (K ^l^L^+,R^l,L^-,R + Vs ^Il^+,R^I,r^-,l) + h.C. (13) 

In the first (second) term in 6H2 both kinds of particles are incident from same (different) 
direction. Thus the spin(charge) momentum is unaffected and hence, as we shall see, the 
corresponding bosonised operator depends only on charge (spin)-fields. 

We will study the effects of both the 2kF and Akp scattering using bosonisation and 
renormalisation group. The fermions may be bosonised using 

i;Rn(x) = hm - — ^e^*^^^«/^ = lim - — ^e'^^^^"^). 



where a is a regularisation parameter, which is finally set to zero. (We ignore Klein factors, 
because we are using the notation of Shankar|T5[| which does not require the use of Klein 



factors. For an explicit comparison of different bosonisation schemes, see Ref. [|16|). Thus, 
the quadratic and the quartic fermion operators are bosonised as 

0± = 4^^±,^ = -^e^v^*± = -^e^-^(^^±^^) (14) 



and O,/, = ^Pl,^lj+,n^l,/n^-,R/L = (^^^^^^^^^^^ (l^) 

Note that as we had expected, the two-fermion operators depend both on charge and spin, 
whereas the four-fermion operators involve either the charge or the spin. 

The renormalisation group flows of these operators to leading order can be found by 
computing their dimensions. Since the operators depend on the flelds, their dimensions 
are conveniently computed using (p as the coordinate field. For the bilinear operators, the 
dimensions are computed as 

<Ot(r. 0)0,(0. 0)>=p^(-^) (16) 

whereas for the operators that are quartic in the Fermi fields, we get 

<Ol,(r,0)O„M>=^J,-^f. (17) 

Clearly, the correlation functions diverge as r~^ for 0± and r""^ for Oc/s as r — > 0. Thus, 0± 
are marginal and Oc/s are irrelevant, just as they are in the noninteracting case. Generically, 
once the barrier becomes stronger, the above leading order analysis of the RG fiow is no 
longer valid, and it is necessary to study the RG fiow to higher orders in V. Since only 
one operator is marginal and the rest are irrelevant, one needs to ask whether the marginal 
operator is marginally relevant or marginally irrelevant when computed to higher orders in 
Ve- But here, it is easy to check that the operator cos 0± can never contribute to itself at any 
order in perturbation theory, because cos"(v^47r0-|-) does not contain a factor cos(v^47r(/)-|-) 
for any n > 1. Hence, the operator in Eq.([T^) or the 2/cf backscattering remains marginal 
to all orders in V^ and all higher order scattering processes {4:kp, Gkp, ...etc.) are irrelevant. 
One may still wonder whether the RG fiow equation for the 2kp backscattering can 
receive non-perturbative contributions. But in this model, since the '+ particles' interact 
only with '— particles' and vice versa, the back scattering of any one kind of particle is just 
not affected by the interaction. Another way of seeing this is to note that this interacting 
system differs from the noninteracting one, by a pseudo unitary transformation which does 
not affect the diagonal density operator ip'lipcr- Hence, even non-perturbative corrections to 
the lowest order flow equations cannot occur. However, we can also see this explicitly if, 
following Kane and Fisher]^, we study the model perturbatively from the other limit - ie., 
the case where there is an infinite barrier with no transmission across the barrier. Then, 
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the fact that the barrier is not infinite, is modelled by a hopping term across the site. We 
then use perturbative RG to determine whether the hopping term is relevant leading to a 
'healing' of the wire, is irrelevant confirming the stability of the infinite barrier fixed point 
or is marginal. 

In first quantised language, for an infinite barrier, the wave function ip^j has to vanish 
at the impurity site x = 0. This allows only odd parity states {\[^„{x) — ipcj{—x)]). In 
terms of the bosonised fields, the condition of infinite barrier at the origin is imposed in the 
following way. We bosonise the < x < L/2 section by introducing the bosonic fields 0q,> 
and 6'a>, and the > x > —L/2 section by the fields 0q,< and 6a<- (The ring boundary 
condition is satisfied by demanding continuity of charge and current at —L/2 = +L/2.) The 
Hamiltonian for the infinite barrier case is then given hj H = J2i Hi, with i =>, < where, 

+ 26{^J.,J29.(|)+:i-9x0+^J29.<t>-J)}■ (18) 

j j 

with the boundary conditions 0<(O) = 0>(O) = | [|1^. Although, we have effectively 
decoupled the wire by allowing for two different kinds of bosons < and > to the left and 
right of the origin respectively, note that the J2j dx(t>aj terms in the Hamiltonian, ensures that 
the '+' particles in any one section (> or <), interacts with all '— ' particles( < and >) in the 
ring|l^. This would certainly not have been true if we had two disconnected wires(rings). 



Since the interaction energy is not extensive, the original ring can never be reproduced 
starting from two completely disconnected rings. The above Hamiltonian, in fact, models a 
single ring with no current across the impurity. (Here we wrote the Hamiltonian in a = ± 
representation, but equivalently, we could have written it in terms of the spin and charge 
fields.) A redefinition of both the < and > fields analogous to the field redefinitions in 
Eq.(^, brings it to the apparently noninteracting form given by 

H=—Y.j dx{{d,c^„,f + {d^e^^f}. (19) 

with a = ± or c, s and i =<, >. 

We now introduce a hopping term across x = 0, as 5H = tJ^aOa, where, 

O^ = ^t (o)^^^(o) + h.c = ^L^e*v^(^->(0)-e<^<(0)) + he 

(27ra) 

= .^_e^V^i9^>W-o.<m) ^ h.c (20) 

(27ra) 

and t is an overlap matrix element. Note that since the hopping operator depends only the 
6cr{0) (since 0o-(O) = 0), we may compute its dimension using the Hamiltonian in Eq.(|l^) by 
treating the 6 fields as the coordinate fields. The dimension of the hopping operators are 
easily computed as 

<Oi(r. 0)0,(0. 0)>=pl^(^), (21) 



which says that the operator is marginal as expected, and hence the single particle backscat- 
tering process (2^^?) does not depend on the interaction. As for the non-interacting case, it 
depends on the strength of the scattering potential and the energy of the incident particles. 
In this limit too, one can show explicitly that the higher order processes are irrelevant. 

Hence, the result that we obtain for transmission when the wire has a non-magnetic im- 
purity at the origin introduced either through a weak 6- function barrier or in the opposite 
limit, through weak hopping, is that the strength of the non-magnetic impurity V is unaf- 
fected by the interaction. The system thus behaves like a non interacting electron system. 
For such impurities, depending on the values of the energy and V, a part of the incident 
electrons gets transmitted and a part gets reflected. Neither perfect transmission nor perfect 
reflection is observed. 

As an aside, note that the dimension of the same hopping operator can also be computed 
in a system of two completely disconnected rings. The Hamiltonian for such a system is 
the same as Eq. (]T^), except that the J2jdx(paj is replaced by dx4>ai- In this model, the 
dimensions are given by 



a' 



<0t(.. 0)0.(0. 0)>=p^(j^j^)«, (22) 

which implies that the operator is irrelevant. So the model with two disconnected rings 
does not get healed; rather it flows to the insulating fixed point where the two rings are 
disconnected. 

What about impurities which can flip the spin of the electron at the impurity site? 
Such magnetic impurities can be introduced through spin dependent potentials. For a weak 
barrier, we may introduce, 

6H = V{^i{0)^40) + h.c). (23) 

At the lowest order, this will generate two different 2kF backscattering processes, iI)\i^iI)-^r 
and ^^^V^__i, both of which can be written in the bosonised form, as, 

O^ = V^l,i^^,R = e^^^('^±+'^T+e±-eT) (24) 

where 0± is just the notation for the operators ipl_ j^i/j^^r and 'ijj_ ^ip+^R respectively. Since 
for computing dimensions, it is more convenient to use the apparently non-interacting form 
of the Hamiltonian given in Eq. (|I0|) , these operators may be rewritten in terms of the tilde 
fields as 

Q^ = gi0F[(l±5)^±+(l±5)0:p-0±+0:p) (25) 

using Eq.(|D The lowest order contribution to the RG equation for V is now simply obtained 
by computing the dimension of the operator. The correlation function is now given by 

<0I(., 0)0,(0, 0)>^(^(^)Jl-<--n (26) 



and diverges as r ^ "^ ^^"^ as r — > 0. So the operator 0+ = V^^ l'^-,r is relevant for (5 < 0, 
whereas the operator 0_ = V^l l'4'+,r is relevant for 5 > 0. Thus one of the two backscattering 
spin flip operators is always relevant for either sign of effective interaction strength. Hence, 
the magnetic impurity cuts the wire at the site of the impurity. Since the lowest order 2kF 
backscattering is relevant, we need not look for higher order processes. An explicit check, 
however, shows that all such operators are irrelevant. Once the impurity strength grows 
under renormalisation, the weak barrier analysis is no longer valid. As was done in the case 
with non-magnetic impurities, we need to study the strong barrier limit, where we start from 
the insulating limit and then study weak hopping. We take the Hamiltonian of Eq. (ITB|) in 



the 6 representation, (or Eq.([T9|) in the 6 representation) which is effectively cut at the site 
of impurity as the unperturbed hamiltonian. The hopping of electrons across the the weak 
link, accompanied by a spin flip is modelled by 

(5if = t(4>(0)^_<(0) + /i.c.), 

which can be bosonised as follows, 

f Of ~ ~ - ~ 

5H = __e^v^(^+>-^-<) + h.c = ^—cos[V^ie^^ - L^ - 6J2k^)] (27) 

f7,2 

where the apparently non-interacting tilde fields have been used in the second line. The 
scaling dimension can easily be read off as 2 + 46'^. So hopping remains irrelevant for all 
values (both +ve and -ve) of the interaction parameter. Here again, we can explicitly check 
that all higher order processes continue to be more irrelevant. The wire which is cut at the 
weak barrier limit, does not get 'healed' for any strenght of the barrier. This interacting 
electron system, thus, flows towards an insulating fixed point in the presence of magnetic 
impurity. 

Let us contrast these results with the results obtained in a purely TR invariant model 
of a Luttinger liquid, - e.g., the model studied by Kane and Fisher^. They found that for 
spinless fermions, repulsive interactions imply that any barrier cuts the wire and attractive 
interactions imply that any barrier is rendered invisible in the renormalisation group sense. 
However, once both spins and charges are introduced with arbitrary values of gp and g^ 
denoting the strengths of the interactions of the charge and spin sectors, various phases do 
exist. Besides the purely conducting and the purely insulating phases, mixed phases where 
there exists a non-trivial fixed point depending on barrier height separating the two phases 
also exist. Kane and Fisher also showed how to access the non-trivial fixed points by studying 
the renormalisation group equations to second order in the impurity strengths around the 
weak and the strong barriers. 

This model cannot be fitted into the general framework of the Kane-Fisher models of 
Luttinger liquids, because it explicitly violates TR symmetry. However, at a qualitative 
level, the results are similar to the Kane-Fisher results [§] for repulsive interactions. At least 
for spin-flip impurities, one finds that any barrier, however small, cuts the wire. For non- 
magnetic or spin-conserving impurities, however, we found that this model is similar to a 
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non-interacting model, and allows both transmission and reflection. Note, however that we 
are not assuming phase coherence through the whole ring as exists for mesoscopic systems. 
In this model, the strengths of the interactions of both the spin and the charge degrees 
of freedom are completely fixed by the gauge interaction. We have no freedom in changing 
the relevant strengths of the interactions, or in fact, explicitly changing the strength of the 
TR violation. However, in a more general model of transport which allows for arbitrary TR 



symmetry violation|]18[, we can try to look for other phases such as the purely conducting 
phase or mixed phases separated by non-trivial fixed points. It should also be possible in 
such models to gradually change the strength of the TR violation and show that in the limit 
where the TR violation goes to zero, the Kane-Fisher results are recovered. Another useful 



generalisation [T^ is to study transport in many-chain (band) models which can be modelled 



for some purposes by a multi-band Schulz-Shastry modelO, 0. Here, again, it may be 



possible to have more phases and perhaps non-trivial fixed points separating the phases. 

In conclusion, we have studied transport in a a simple exactly solvable model of in- 
teracting fermions and have obtained the anomalous dimensions of several back-scattering 
operators in the weak barrier and weak link limits. Since at low temperatures and low driving 
forces, the conductance has a power law behaviour in terms of the anomalous dimensions, we 
can directly obtain results for transmission and conductance. With several new experiments 
on one dimensional wires (semi-conductor wires in the single channel limit, single wall carbon 
nano-tubes, etc.), it becomes essential to include interactions. The usual Landauer-Buttiker 
formalism will not suffice. Hence, it is important to study transport in models with inter- 
actions. Besides the results that one gets in the standard Luttinger model, it will definitely 
be extremely useful to have other models of interacting fermions, where many of these ideas 
of transport through LLs in various geometries can be tested. Since any external voltage 
applied to the wires breaks TR invariance, it is also of interest to study transport in models 
where TR symmetry is explicitly violated. In this context, we feel that the exactly solvable 
models proposed by Schulz and Shastry could serve as a useful tool. 
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